from gruppenFinding import * 
	
	
cutoff = 1
proteinMapping = readDictionary()
conservations = buildConservations(cutoff, proteinMapping)
inpalist = {}
species = "Bacteroides"

for protein in conservations:
        if str(protein) in proteinMapping:
                if proteinMapping[str(protein)]==species:
                       # print("prot von bacteroides")
                        if species in protein.specMatches:
                                #print("bt in specmatches")
                                inparaloge = [str(match) for match in protein.specMatches[species]]
                                inparaloge.append(str(protein))
                                inparaloge.sort()
                                inpalist[str(inparaloge[0])] = str(inparaloge)

print(len(inpalist))
                                

                                
			
